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ABSTRACT 



Different ray paths through a turbulent plasma can produce stochastic Faraday 
rotation leading to depolarization of any linearly polarized component. Simple 
theory predicts that the average values of the Stokes parameters decay according 
to (Q), {U) oc exp(— (5;), with 5i cx A^. It is pointed out that a definitive test 
for such depolarization is provided by the fact that (Q^ + C/^) remains constant 
while (Q)^ + (C/)^ decreases oc exp(— 25;). The averages to which this effect, called 
polarization covariance, should apply are discussed; it should apply to spatial averages 
over a polarization map or temporal averages over a data set, but not to beamwidth 
and bandwidth averages that are intrinsic to the observation process. Observations of 
depolarization would provide statistical information on fluctuations in the turbulent 
plasma along the line of sight, specifically, the variance of the rotation measure. Other 
effects that can also cause depolarization are discussed. Favorable data sets to which 
the test could be applied are spatial averages over polarization maps of the lobes of 
radiogalaxies and temporal averages for the binary pulsar PSR B1259-63 as it goes 
into and comes out of eclipse by the wind of its companion. 
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1. Introduction 

Faraday rotation can cause depolarization of an astrophysical synchrotron source in three 
ways: internal depolarization, bandwidth depolarization and external depolarization (e.g.. Burn 
1966). Internal depolarization is due to differential Faraday rotation when a ray emitted at the 
far side of the source has its plane of polarization rotated through ^ vr/2 before reaching the near 
side of the source. Bandwidth depolarization occurs due to differential rotation of the plane of 
polarization across the bandwidth of observation (e.g., Simonetti, Cordes & Spangler 1984; Lazio, 
Spangler & Cordes 1990). External depolarization is due to a random component of the magnetic 
field along the ray path between the source and the observer through the intervening medium 
(interplanetary, interstellar, intergalactic or intracluster) . Two different treatments of external 
depolarization are available. One treatment, referred to here as stochastic Faraday rotation, is 
based on the assumption that the (Faraday) angle through which the plane of polarization is 
rotated contains a random component, (e.g., Burn 1966; Simonetti et al 1984; Tribble 1991). 
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The other treatment of external depolarization involves including the effect of the magnetic field 
in scintillation theory (Erukhimov & Kirsh 1973; Tamoikin &: Zamck 1974; Melrose 1993a, b 
[hereinafter papers A and B]). Simonetti et al (1984) established that these two treatments arc 
equivalent for the mean values of the Stokes parameters, as shown in the Appendix below. Such 
stochastic Faraday rotation is of practical interest in determining the properties of the medium 
through which the radiation propagates: measurement of the frequency at which depolarization 
occurs provides information on the fluctuations in the magnetic field along the line of sight 
which complements information from the dispersion measure (DM), rotation measure (RM) and 
scattering measure (SM), e.g., Taylor &: Cordes (1993). However, in cases where depolarization is 
observed it is not necessarily clear that it is due to stochastic Faraday rotation, rather than to one 
of the other depolarizing mechanisms. 

In the present paper, it is shown that a definitive test for stochastic Faraday rotation follows 
by measuring the mean square values of the Stokes parameters, specifically (Q^), (C^^), {QU). 
These mean values can be derived in a very simple way by assuming that the Faraday rotation 
angle can be treated as a random variable. Let the Faraday angle be 

= <Po + RMX\ (1) 

where 0o is the angle at the source and A is the wavelength. The statistical average is performed 
by assuming that (p-p has a random component that obeys gaussian statistics. This simple model 
is justified in an Appendix, where it is shown that the theory of scintillations (e.g. Narayan 
1992) in a magnetized plasma leads to the same result for (Q^), {U"^), {QU) as this model. The 
model implies that (Q^ + U'^) is a constant, that is, although (Q), (U) decay, (Q^ + U^) does not. 
This simple but surprising result is referred to here as polarization covariance; it distinguishes 
depolarization due to stochastic Faraday rotation from other depolarization mechanisms. 

In section 2 a model for stochastic Faraday rotation is used to derive formulae that describe 
the depolarization and the polarization covariance. In section 3, the nature of the averaging 
process is discussed. In section 4, some specific examples of possible depolarization are considered. 
The conclusions are summarized in section 5. 

2. Random Variations in the Rotation Measure 

Faraday rotation results from the difference between the refractive indices of the two natural 
modes of a magnetized plasma. In the simplest approximation the polarization of the natural 
modes is assumed circular, and then Faraday rotation affects only the linear polarization. Here the 
polarization is described in terms of the Stokes parameters, /, Q, U, V. Only Q, U are affected 
by Faraday rotation; I, V are not mentioned further in this paper. 
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2.1. Faraday rotation 

Faraday rotation causes the Stokes parameters Q, U to vary with distance z along the ray 
path according to 



dz \U{z) J ' ' Vl Oj \U{z)J ' 

It is convenient to write equation (^) in the form (e.g., Spangler 1982) 

^V{z)=ipAz)V{z), V = Q + iU. (4) 
The solution of equation (Q) is 

Mz) = r dz' py{z'), (5) 



where 7^0 = ^0 + ^Uq is determined by the intrinsic polarization at the source (at z = 0). In the 
following it is assumed that the axes are oriented such that Uq = 0, so that one has Vo = Qo- 



2.2. Fluctuations in the Faraday angle 

Fluctuations in the plasma parameters along the ray path may be taken into account by 
separating the Faraday angle, (pp, into an average value, (</>f), and a fluctuating part, 5(f)-F- 

(pF = (0f) + ('?^'f) = {Pv)L, (6) 

where L is the length of the ray path. It is convenient to incorporate the average Faraday rotation 
into modified Stokes parameters, denoted by tildes (paper A): 

V{z) = e-'^^^^V{z) = e'^^^'^'^^Qo, (7) 

where the final equality follows from equation (|5|). One has V = Q + iU, with Q, U real. The 
fluctuating part of the Faraday angle, ScpY, on the right hand side of equation ^ is now assumed 
to be a random variable. 



2.3. The mean Stokes parameters 



The mean values of the Stokes parameters follow from the statistical average of equation 
d?]), which may be evaluated using the property, (exp(i(/))) = exp(— ((/)^)/2) for a gaussian random 
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variable, (/>. (The derivation of this property follows from the power series expansion of the 
exponential and {<f>^'^'^~^^^) = 0, (0^") = 3 • 5 • . . . • (2n — for n > and integer.) This gives 

{V) = QoeM-kii^^Ff)]- (8) 
It follows that the Stokes parameters vary according to 

(Q) = Qoe-'S (^7) = 0, (9) 

where equation (Q) is used and 6i = ((5(/>_p)^)/2. The decay of {Q) described by equation (P) 
corresponds to a depolarization of the radiation. This result has been derived both by the foregoing 
method (Burn 1966, Tribble 1991) and also from scintillation theory for a magnetized plasma in 
the special case of axial rays (Erukhimov & Kirsh 1973; Tamoikin & Zamek 1974; Simonetti et al 
1984; paper A). The exponent of the decay factor in equation (^) is determined by 

((5#)2) = r dz ("^ dz' {5py{z) 8py{z')), (10) 



with 5pv = Pv — {Pv)- Equation (10) may be rewritten by analogy with equation (^) in the form 
((5<^f)2) = {Spl)L^. 



2.4. The mean square Stokes parameters 

The mean values of the squares of the Stokes parameters and of the products of the Stokes 
parameters follow from the statistical average of the modulus squared of (^) and of the product 
of equation (|^) and its complex conjugate. Under the assumption that dcpp is a, gaussian random 
variable, these give 

(^2) = Qlexp{-45i), 
(Pv*) = Ql (11) 

In terms of the Stokes parameters, equation (|l|) implies 

(Q') = cosh(250, 
(U^) = Qle-^^'smh{25i), 

(QU) = 0. (12) 



The result (12) was derived in paper B by solving the fourth order moment equation for a 
magnetized plasma for axial rays. The fact that the foregoing theory reproduces the result of the 
more general theory confirms the equivalence of the two theories for the higher order moments for 
axial rays. 
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The covariances of the polarized component of the radiation are given by 



cosh(2(5;), 

sinh(2(5z). (13) 



It follows that for small 5i the mean square of Q is equal to the square of its mean; however, as 5i 
increases, the ratio of the mean square fluctuations to the square of the mean increases, becoming 
arbitrarily large for 5i ^ 1. The mean of U is zero by definition, and the mean square of U is 
small compared to the mean square of Q for small 5i, but one has {Q^) ~ {U"^) for 5; ^ 1. This 



somewhat surprising prediction may be attributed to the second of equations (lll|) : because the 
Faraday angle does not appear, the value of + U"^ is independent of the statistical averaging. 
Whether or not this is the case for averaging processes in practice is discussed in the next section. 



2.5. Consequences of polarization covariance 



The relations (|8|), (|11|) and may be rewritten in a variety of ways. The tildes on the Stokes 
parameters correspond to, cf. equation (|^) 

Q = Qcos(0f) + f^sin((^F), 

il = -Qsin(0F) + t/cos(0F). (14) 

The relations (12) involve the parameters (0f) and 5i, but not Qq. Various relations between the 
averages of the Stokes parameters (without the tildes) may be derived by using equation ( [l^ ) and 
eliminating one or more of these three parameters. For example, relations that involve Qo) but 
not ((/>f), are 

m + {U^) - (Q)2 - = Ql[l - e-'''), (15) 

{Q^){U')-{QUf = lQl(l-e-''^). (16) 

More generally, there are five observable parameters in this theory, {Q), {U), (Q^), {W^), {QU), 
and three unknowns, {(fip), Si, Qq. Thus measurement of all five observables would provide two 
constraints on the theory. One of these is equivalent to (U) = implying (QU) = 0, and the other 
is a definitive test for stochastic Faraday rotation. 



3. Conditions for Polarization Covariance to be Observed 

The conditions under which polarization covariance should be observed may be identified 
by considering the manner in which the observational data are recorded and the nature of the 
averaging process intrinsic in, or applied to, the data. For a given radiotelescope observing at a 
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given frequency, data are recorded over an integration time, At, and over a number of frequency 
channels, each of bandwidth, Ai^, centered on the given frequency. In principle, the basic bits of 
information are the values of /, Q, U (and maybe V) for each At and Ai^. Over an observation 
time, the averages of /, Q, U are simply the averages of the values for each At, and similarly the 
values of I^, Q^, C/^, QU, etc., are to be obtained by storing the appropriate products for each At 
and then averaging each of them over the observation time. The data may also be averaged in 
frequency over all the channels. Furthermore, the telescope has a characteristic beamwidth, A9, 
and the bits of data arc already averaged over this beamwidth. Hence, one may identify three 
types of average: intrinsic averages in the data over At, Av, AO; the averages applied to the 
data over the observation time and bandwidth; and temporal averaging over times long compared 
with a given observation or spatial averages on scales large compared with the resolution of the 
telescope. 

3.1. Averages over space or time 

The interpretation of an ensemble average, denoted by the angular brackets above, is not 
specified a priori. It should apply to the averages over space, time, or beamwidth. As the space 
scale or time scale over which the average is performed is increased, the value of 5i should increase 
(or remain constant). In principle, measurement of 5i provides information on the fluctuations 
over the relevant time or space scales. 

Spatial averaging requires a map showing the variation of RM (e.g., Simonetti et al 1984, 
Simonetti &: Cordes 1988). To simplify the discussion, suppose that the Faraday rotation all occurs 
in a screen at a distance L. An angular separation, 89, corresponds to a transverse displacement 
r = L60 at the screen, and hence an average over 59 corresponds to an average over transverse 
displacement at the screen. Simonetti et al (1984) introduced a structure function for RM: 

UKum = {[RM{9) - RM{e + 69)f), (17) 

and showed that it consists of a statistical part and a geometric part. The statistical part for a 
power law spectrum of fluctuations is of the form 

KKMm oc {59f (18) 

for < /3 < 2, with /3 = 5/3 corresponding to a Kolmogorov spectrum. The geometric part is due 
to the different path lengths through the random medium, and it appears to be unimportant in 
practice. 

The condition for the theory of stochastic Faraday rotation to apply is that the variations 
in the phases of the components in the two modes be large, so that their difference, and hence 
the Faraday angle, may be treated as a random variable (Lee & Jokipii 1975; Simonetti et al 
1984). The phase difference between different points would be known for an ideal map, but in 
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practice the sampling is coarse and the separation over which phase coherence is lost is not known 
a priori. The theory for depolarization, as described by equations (§) and (|l^), applies only over 
separations large compared with that over which phase coherence is lost. 

Temporal variations in the Faraday angle can arise from two effects: intrinsic temporal 
changes in the screen (e.g., due to wave motions therein), and convective motions, involving 
relative motions between the the screen and the observer. Only convective motions are considered 
here. An angular displacement, A^, of the image then occurs in a time interval. At, related by 

LA9 = vAt, (19) 

as the pattern sweeps across the observer at speed v. 

For an unresolved source with an apparent size O^s, the minimum time on which different ray 
paths are sampled is 

tmin = min [0cs, Ob]L/v, (20) 

and temporal averaging is relevant only over times ^ tmin- For a resolved source that fills the beam 
the minimum time for temporal averaging is set by the beamwidth of the telescope, imin = L6b/v. 
By way of illustration, for a screen at a distance of several kiloparsecs and with v ~ 100 km s^^ 
(of order the velocity of the Earth through the Galaxy), even for a source of apparent angular 
size ^eff ~ 1 mas, equation ( [20| ) implies that imin ~ 2 months is the minimum time for a temporal 
average to show the effect of stochastic Faraday rotation for observations with milliarcsecond 
resolution. The time scale for temporal variations is shorter for observations with higher angular 
resolution and for closer screens. The effect is of no practical interest for extragalactic screens. 



3.2. Beamwidth averaging 

Beamwidth averaging is relevant for an extended, resolved source, that is, for a source with 
angular size, 9s, satisfying 6s > Ob- A source with 9s <^ Ob is unresolved and beamwidth 
averaging is not relevant. 

The observed intensity for a resolved source is the actual intensity convolved with a point 
spread function, p{r), where r is the perpendicular displacement on an image screen (e.g., Rohlfs 
& Wilson 1996, Eqn. (5.64)). For a gaussian beam one has 

p{r) = exp(-rV2ro2), (21) 

where rg is a constant. 

Consider an idealized model in which S(f)-Fi^,z) varies due to a wave, with wavenumber K 
and amplitude 5B, propagating along the x-axis with magnetic fluctuations along the z-axis. 
Specifically, assume 5(j)F{r,z) = fosinKx, with /o = {6B / Bii)RMX^ . If the beamwidth were zero 
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one would observe V = Qo exp(i/o sinETx). On convolving this with the point spread function 



(21), one finds 

oo 

{V)b = Qi^ E Jn(/o)e*"''"e-"'^''-o/2, (22) 

n=— oo 

where the subscript b denotes the average over beamwidth. In the limit tq ^ the final 
exponential in equation (p^) goes to unity, and the sum gives V = Qo exp(z/o siniTx), as required. 
For Kr^ ^ 1, {V)h decreases sharply and only the term n = contributes to the sum in equation 

m 



Although this model is highly idealized, it suffices to illustrate that significant depolarization 
occurs for Kvq ^ 1, that is, when the fluctuations in RM occur over a scale, ~ 1/K <^ tq, that is 
small compared with the resolution of the telescope. 

For an unresolved source, 9$ < Ob, the beamwidth average is replaced by an average over the 
apparent source. The apparent size cannot be smaller than the scatter-image size of a point source 
affected by scattering in the interstellar medium (ISM), or other turbulent medium through which 
the radiation passes (e.g., Rickett 1990). A simple model for the angular broadening due to such 
scattering (e.g., Blandford & Narayan 1985) involves a gaussian profile similar to that assumed 
above in performing the beamwidth average. Now /(r) oc exp(— r2/2ri^) is interpreted as the 
brightness distribution as a function of r for an apparent source of angular size ~ ri/L. Thus, for 
a point source the average over beamwidth is effectively that over the angular size of the scatter 
image, rather than over the beamwidth of the radiotelescope. 



3.3. Average over bandwidth 

An average over bandwidth causes depolarization because (pp oc u^'^ is a function of frequency. 
Over Ai' this causes a change |A0f| = |2(/)fAi//i/|, and hence a fractional reduction in the degree 
of linear polarization by 2Au/u. Depolarization due to bandwidth averaging is qualitatively 
different from the other averages considered above because it affects the amplitude rather than 
the intensity. Polarization covariance applies to the statistics of the measurements of Q and U, 
and bandwidth averaging limits the ability to measure values of Q and U. 



4. Linear depolarization measure (LDM) 

Assuming that depolarization is observed due to stochastic Faraday rotation in the ISM, the 
wavelength at which the depolarization becomes important provides information on the statistical 
properties of the ISM. This may be described in terms of a 'linear depolarization measure', LDM, 
defined, cf. equation @, 

Si = i((<5(/)F)') = |LDMA^ 
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LDM = ((RM- (RM))2), (23) 

where equation (Q) is used. The quantity LDM (hnear depolarization measure) is determined by 
the variance in RM, and may be measured by any observation of stochastic depolarization. 

Four observationally determined parameters have been used to describe the statistical 
properties of the ISM: the dispersion measure, DM, the emission measure EM, the scattering 
measure SM and the rotation measure RM (e.g., Cordes & Lazio 1991; Cordes et al 1991). These 
parameters determine the integral along the line of sight of, respectively, Ue, n^, Cfj and nei?||, 
where Cfj is the structure constant for the electron density fluctuations. The additional parameter, 
LDM, is determined by the integral along the line of sight of (ree-By — (nei?||))^. 

To illustrate how the LDM could be used to infer information concerning fluctuations in the 
magnetic field, consider two extreme cases (e.g., Simonetti & Cordes 1988). On the one hand, if 
the magnetic field in the ISM is relatively uniform, then the fluctuations in (ne-By)^ are dominated 
by those in ng, and one would expect the ratio LDM/RM^ to vary across the Galaxy in the same 
way as does the ratio SM/DM^. On the other hand, if the magnetic field is sufficiently random, 
or if regions of opposite sign of B\\ nearly cancel, then one can have (-By) ^ (-^||)^) implying that 
LDM/RM^ could exceed unity, due to near cancellation of contributions to RM of opposite sign 
along the line of sight. 

Existing data suggest that the variance in RM in the ISM is of the same order but not much 
larger than (RM)^ (Simonetti et al 1984; Simonetti & Cordes 1988). The implication is that i?| 
does not change sign many times along these lines of sight. Further data confirming (or otherwise) 
this result is clearly important for our understanding of the structure of the magnetic field in 
the ISM. In principle, polarization covariance can be used to estimate (RM)^ with much finer 
resolution than can be achieved by averaging over a polarization map. 

5. Observations of depolarization 

Depolarization observed in some sources has been interpreted in terms of stochastic Faraday 
rotation, but the definitive test for polarization covariance has not been performed. In this section, 
some specific examples are discussed. 

5.1. Depolarization in the ISM 

Lazio et al (1990) used polarization data for eight double-lobed radio galaxies to discuss 
variations in RM in the Cygnus region. They reported significant depolarization, for several of 
these sources, at A ^ 10 cm. For sources for which RM could be determined, values in the range 
~ 250-800 rad m~^ were found. The data were compared with a model that gives a relation of the 
form, cf. equation (0), K{e) = Koie/Bof, with kq ~ (10radm-2)2 for 0q ~ 10" and with p = 5/3. 
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The data are described moderately well by the model, but a smaller value of /3 ~ 1 would appear 
to give a better fit. 

The wavelength, Aq, at which depolarization due to stochastic Faraday rotation should occur 
follows from equations (§) and @: 

.1/2 \ X-/3/2 

where 6b is the beamwidth of the telescope. With Kq "^ = lOradm"^ for a beamwidth of ~ 10", 



.1/2 



equation (24) implies depolarization at A 30 cm. Thus, stochastic Faraday rotation is plausible 
for the depolarization reported by Lazio et al (1990). Such depolarization would be due to 
variations in RM on a scale r ^ L9b, below the resolution of the telescope. 



5.2. Depolarization in intergalactic media 

It has been suggested that differences in the degree of polarization in the two lobes of radio 
galaxies might be due to either internal depolarization or due to stochastic Faraday rotation in a 
circumgalactic medium (Strom & Jagers 1988; Laing 1988; Garrington, Leahy & Conway 1988). 
The observed depolarization occurs typically at A between 6 cm and 20 cm. An example of an 
extragalactic source for which detailed polarization data are available is the inner 2 kpc of M87 
(Owen, Eilek & Keel 1990). This source has relatively high polarization, ^ 10%, in the range 
4.6-4.9 GHz, with large values of RM ~ 1000-8000 radm'^ that vary over scales ^ 1 kpc across 
the source. 

Suppose that there are fluctuations in RM on a scale smaller than can be observed, with 
variance LDM. Then the frequency at which the depolarization should become significant is 

^ lOGHz [LDM^/^/1000radm~^]. The absence of substantial depolarization indicates that 
LDM^/^ cannot be much larger than lOOOradm"^. However, it would be surprising if the magnetic 
field were so uniform that one had LDM^/^ <C RM, and hence one should expect depolarization 
to become important at frequencies of order of those already observed. Tribble (1991) used the 
theory for beamwidth depolarization to place a limit on the value of B in the circumgalactic 
medium. 

The data for M87 should allow one to perform spatial averages of Q, U and + U'^ over 
maps at the available frequencies, and this would provide a test for the interpretation of the 
depolarization, which should affect Q, U but not + C/^. However, this test has yet to be made. 
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5.3. The Snake 

The Snake (Gray et al 1991) is a long (~ 20'), thin 10") structure roughly perpendicular to 
the galactic plane at ~ 1° from the Galactic Centre. It is highly polarized, ~ 60%, at 10.55 GHz, 
and less polarized at 4.79 GHz (Nicholls & Gray 1992). There is an abrupt decrease in polarization 
roughly at its mid point (Gray 1993). 

For the decrease in polarization to be due to stochastic Faraday rotation requires an abrupt 
increase in the fluctuations in RM at the point where the polarization decreases. One possibility is 
that there is an interstellar cloud along the line of sight covering only the weakly polarized portion 
of the Snake. Testing this suggestion requires multifrequency observations to determine RM along 
the Snake. 

5.4. Pulsars 

Most pulsars are highly polarized (e.g., Lyne & Manchester 1988), implying that stochastic 
Faraday rotation is normally unimportant. However, there are exceptions which could provide 
useful tests for the theory of stochastic Faraday rotation. 

First consider the frequency at which depolarization might be expected due to variations 
in RM in the ISM. The relevant variations are over the apparent angular size, which along 
with the time delay and the pulse width may be derived from a simple model for refractive 
scattering (Blandford & Narayan 1985, Rickett 1990). Gwinn, Bartel & Cordes (1993), found 
that the angular size is reasonably well fit by an estimate based on the temporal broadening, r, 
through Or = 3.3(ct/L)^/2. At 327 MHz angular sizes ranged from the detection limit(~ 2 mas) 
to ~ 65 mas. Using equation (p^ with 9b replaced by the apparent angular size, ^ 100 mas, one 
concludes that depolarization for most pulsars should occur only at unobservably low frequencies, 
consistent with the observations. 

One exceptional case where depolarization of pulsars might be observable in a region around 
a supernova remnant (SNR) where Alfven waves are excited in association with shock acceleration 
of relativistic particles (Spangler et al 1986). Gwinn et al (1993) argued that such turbulence may 
account for strong scattering around young SNR. These authors noted that there are variations 
in DM with time up to 0.01-0.06 pccm~^. Assuming a magnetic field ~ 3/^G, this suggests 
the contribution to RM from the turbulent region around the SNR could be comparable to the 
total RM ~ 100-1000 radm^^. In a model for SNR with sharp edges, Achterberg, Blandford Sz 
Reynolds (1994) suggested that the required Alfven waves imply fluctuations in the magnetic field 
SB/B ~ 3 X 10~^ on a scale ~ 10^*^ m, which is smaller than the scattering disc. This model 
implies fluctuations in RM across the scattering disc as large as ~ Iradm"^, which suggests that 
depolarization should occur at A 1 m. 

A clear example where depolarization does occur is for the binary pulsar PSR B1259-63 which 
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is eclipsed by the wind of its companion for several weeks near periastron. Large fluctuations 

in DM and RM were observed in the 1994 eclipse (Johnston et al 1996), and, in the following 
eclipse in 1996, depolarization was observed over a few days (around day —80) before the pulsar 
disappeared and for several days as it came out of eclipse (S. Johnston, private communication 
1998). This depolarization was clearly associated with fluctuations in RM on a time scale of 
minutes. This system provides the possibility of testing the theory of stochastic Faraday rotation 
in detail. Multifrequency observations of the polarization on a short time scale in the epochs 
where the variations in RM occur would provide a data set allowing the averages oi Q, U and 
q2 ^ jj2 performed over various time scales as the level of fluctuations in RM changes. It is 

hoped that such a data set will be recorded around the next periastron. 

6. Conclusions 

It has long been known that stochastic Faraday rotation can lead to depolarization (Burn 
1966). Writing Q = Qqcos^f, U = Qosin^p, for the observed Stokes parameters in terms of 
Qo at the source, and assuming <pp = (f)Fo + (50f with S(f)F obeying gaussian statistics, implies 
(Q) = QoCOS(/>Fo exp[— (5j], (^7) = QosincppQ cxp[— (5/]. It is pointed out in the present paper that 
this depolarization does not affect {Q'^ + U'^), which should remain constant while (Q)^ + {U)'^ 
decreases. This effect, which is referred to as polarization covariance, was first identified within 
the framework of a more general theory for scintillations in an anisotropic medium (paper B), and 
stochastic Faraday rotation provides an obvious interpretation of it. The interpretation is that 
depolarization due to stochastic Faraday rotation is due simply to random rotation of the Stokes 
vector. In principle, whether depolarization is due to stochastic Faraday rotation or to some other 
effect may be tested by recording data on Q and U on the shortest time and space scales available, 
and calculating + C/^ and (j)p from these data. In stochastic Faraday rotation, (pp varies 
randomly while + U'^ remains constant. The variance in the rotation measure is proportional 
to {4>p), but the same information is provided by comparison of (Q)^ + (U^) and {Q^ + U^). 
The latter quantities could be recorded automatically on the shortest time scale available on a 
telescope by modifying the software in polarization measurements to record the value of + U^, 
as well as the values of Q and U, so that the integrating over the time during the observation gives 
(g2 + U^)^ (Q) and {U) directly. 

In principle, polarization covariance provides a definitive test for depolarization due to 
stochastic Faraday rotation, and the conditions under which this test should apply are discussed 
in section 3. There arc at least five relevant averages: 

1. Beamwidth averaging: the average is over the beamwidth of the observation. If the source 
is larger than the scattering disc but smaller than the telescope beam, this average is over the 
unresolved source. 

2. Averaging of a scattering-broadened image: the average is over the scattering disc. 
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3. Bandwidth averaging: the average is with respect to frequency over the bandwidth of the 
observation. 

4. Spatial averaging: the average is with respect to position over a polarization map. 

5. Temporal averaging: the average is over a temporal sequence of observations. 

Of these five, the first three are intrinsic to the collection of the data, and polarization 
covariance does not apply to them. However, for spatial and temporal averages, stochastic Faraday 
rotation leads to a depolarization in which [Q)'^ + {U)"^ decreases and [Q"^ + C/^) remains constant. 

Any observation of such depolarization provides information on the statistical properties of 
fluctuations in the magnetized plasma. Specifically, the wavelength at which the depolarization 



sets in may be used to define a linear depolarization measure, LDM, cf. equation (23), which 
provides a measure of the variance in RM, and hence on the fluctuations in (ne-By)^ along the 
ray path. Available data on the variation in RM in the ISM suggest that is slightly 

greater than unity (Simonetti et al 1984; Simonetti &: Cordes 1988). 

The definitive test for depolarization due to stochastic Faraday rotation has yet to be applied 
to any specific data set. Possible systems for which depolarization is plausibly due to stochastic 
Faraday rotation are discussed in section 5. Polarization maps of the lobes of radiogalaxies could 
be used to test the theory for spatial averaging. The next eclipse of the binary pulsar PSR 
B1259-63 should provide an almost ideal test of the theory for temporal averaging. 



Appendix 

The theory of strong scintillations (e.g., Ishimaru 1978) may be developed in terms of 
equations for the moments of the wave amplitude. The wave amplitude, n(z,r), is described 
in the parabolic approximation; with the paraxial direction along the z axis and with r the 
two-dimensional vector perpendicular to this direction. The equation for the nth order moment 
involves a first derivative with respect to z and second derivatives with respect the n rs. The 
statistical average over the fiuctuations is based on an ensemble average assuming gaussian 
statistics. The generalization to an anisotropic medium was made by Kukushkin &; Ol'yak (1990, 
1991) and in papers A and B). In a simple plasma in which the modes are circularly polarized, the 
dielectric tensor, Kij, may be separated into an anisotropic part, Kj, and a part, Ky, associated 
with the anisotropy. Three correlation functions are needed to describe the fluctuations, SKj, 
5Kv- these are written AxY{z,r) = AYxiz,r), which are proportional to {6Kx{z,r)6KY{z,0)), 
with X,Y = I, V. Alternatively, the fluctuations may be described in terms of the phase structure 
function, D{z,r) = {[(f){z,r) - (f){z,0)]'^), with ^^(2:, r) = 2Az[A{z,0) - A{z,r)] for a 'thin screen' 
of thickness Az. The generalization to an anisotropic medium involves three such functions, 
DxY{z,r), constructed from the three AxYiz,r). When considering the effects of Faraday rotation 
on the correlation functions of Q and U one is interested in the case where all the rs are zero. In 
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paper A it was shown that the nontrivial second order moments then give 

,2 



d o;^ , fl 0\f{Q) 

l(^)j+^^^^U i)[{u) 



0, 



Integration of (A.l) over a slab of thickness Az reproduces equation with 

,2 



UJ 



26i = Aw Az. 



(Al) 
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The fourth order moments include the correlation functions between the Stokes parameters. In 
paper B it was shown that, when all the rs are zero, the nontrivial fourth order moment equations 
then reduce to 
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im)\ 

{UU) = 0, 

m)) 



(A.3) 



where the tildes have the same meaning as in section 2, and with Ayy = Ayviz, 0). Integration of 
(^.1) over a slab of thickness Az reproduces equation (12) with 5i given by {A.2). This formally 
justifies the treatment of the Faraday angle (j)^ as a random variable in section 2 in discussing the 
correlation functions between the Stokes parameters. 
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